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ABSTRACT 



Context. In the context of star and planet formation, understanding the formation of disks is of fundamental importance. 
Aims. Previous studies found that the magnetic field has a very strong impact on the collapse of a prestellar cloud, 
by possibly suppressing the formation of a disk even for relatively modest values of the magnetic intensity. Since 
observations infer that cores have a substantial level of magnetization, this raises the question of how disks form. 
However, most studies have been restricted to the case in which the initial angle, a, between the magnetic field and 
the rotation axis equals 0° . Here we explore and analyse the infiuence of non aligned configurations on disk formation. 
Methods. We perform 3D ideal MHD, AMR numerical simulations for various values of fi, the ratio of the mass-to-flux 
to the critical mass-to-flux, and various values of a. 

Results. We find that disks form more easily as a increases from to 90° . We propose that as the magnetized pseudo- 
disks become thicker with increasing a, the magnetic braking efficiency is lowered. We also find that even small values 
of Q (~ 10-20°) show significant difi^erences with the aligned case. 

Conclusions. Within the framework of ideal MHD, and for our choice of initial conditions, centrifugally supported disks 
cannot form for values of fi smaller than ~3 when the magnetic field and the rotation axis are perpendicular, and 
smaller than about ~5-10 when they are perfectly aligned. 
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1. Introduction 

Understanding the formation of protostellar disks is an 
important issue in the context of star and planet forma- 
tion. The presence of circumstellar disks is well established 
around T-Tauri stars (Watson et al. 2007), as well as around 
younger class I and II protostars. However, because of dif- 
ficulties in disentangling the emission of the disk and enve- 
lope, the existence of disks around class objects remains 
a matter of debate (Mundy et al. 2000, Belloche et al. 2002, 
Jorgensen et al. 2007). 

Starting with a level of rotation consistent with obser- 
vations, hydrodynamical numerical simulations discovered 
that centrifugally supported disks form with characteristic 
diameters of the order of a few hundred AU (Matsumoto 
& Hanawa 2003, Hennebelle et al. 2004, Goodwin et al. 
2007, Commergon et al. 2008). These structures are the 
unavoidable consequence of angular momentum conserva- 
tion during the collapse phase. In addition, since these disks 
are massive and strongly self-gravitating, they are prone to 
gravitational instability and have often been found to frag- 
ment, leading to the formation of a small cluster of stars. 

However, magnetic fields appear to have an important 
impact. Magnetohydrodynamic (MHD) simulations (e.g.; 
Machida et al. 2005, Banerjee & Pudritz 2006) find that 
even for modest values of the magnetic intensity, disk for- 
mation can be suppressed by magnetic braking (Shu ct al. 
1987, Mouschovias 1991, Basu & Mouschovias 1995, GaUi et 
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al. 2006) which transports angular momentum from the in- 
ner parts of the cloud towards its outer regions (Allen ct al. 
2003, Fromang et al. 2006, Price & Bate 2007, Hennebelle & 
Fromang 2008 (hereafter HF08), Mellon & Li 2008, 2009). 

The typical values of ^, the ratio of the mass-to-flux 
to the critical mass-to-flux ratio, at which the formation of 
a rotationally supported disk is suppressed, vary from one 
study to another; probably because of different choices of 
initial conditions. For example, while Mellon & Li (2008) 
estimate that the value of /.t, below which disks do not form, 
is larger than 10, Price & Bate (2007) and HF08 find that 
this occurs when /i < 5 — 10. A broad distribution of /i has 
been inferred from observations, but most cores have mag- 
netic fields corresponding to values of fi typically smaller 
than 5 (Crutcher 1999). This raises a fundamental question 
about how disks can form in spite of the relatively high de- 
gree of magnetization. 

Most theoretical studies have been performed by assum- 
ing a simple initial configuration in which the magnetic field 
and the rotation axis are parallel. Notable exceptions are 
the works of Machida et al. (2006), who explored the influ- 
ence of a non zero angle between the magnetic field axis and 
the rotation axis, and Price & Bate (2007) who considered 
the case where the two axes are perpendicular. However, 
none of these works focused on the problem of disk forma- 
tion and magnetic braking in great details, which arc the 
topics of the present work. 

The paper is organised as follows: in sect. 2, we analyse 
the magnetic braking emphasizing the importance of the 
initial angle, a; in sect. 3 numerical simulations of coUaps- 
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Fig. 1. Density and velocity profiles. First column is for /.t = 16 and a = 0°, second column is for \x = h and a = 0°, 
while third column is for /i = 5 and a = 20°. The first row shows the x-y plane (therefore perpendicular to the magnetic 
field direction) and a zoomed image of the centrifugally supported disk (when it exists) . The second row shows both the 
y-z plane and the more extended magnetized pseudo-disk seen edge on (| y |> 200 AU). 



ing cores for various values of a and [i are presented; sect. 4 
concludes the paper. 

2. Magnetic braking: physical analysis 

The classical analysis of magnetic braking, presented in 
Appendix, assumes that there is a clear distinction between 
the cloud and the intcrcloud medium. However, inside a col- 
lapsing core the situation differs a priori. Density gradients 
(typically p oc r^^) throughout the cloud make it difficult 
to discriminate between its dense inner parts and the dif- 
fuse external medium. It is therefore not entirely clear to 
what extent the classical analysis may be applied to under- 
stand the magnetic braking operating inside a collapsing 
dense core. Indeed various studies have already demon- 
strated that internal magnetic braking is more complex 
(e.g.; GaUi et al. 2006, HF08). 

A simple analysis of the terms responsible for angular 
momentum transport through the cloud is elucidating. We 
consider a collapsing cloud permeated by a uniform mag- 
netic field, Bz- Rotation with respect to the z-axis generates 
a toroidal field, Bq (r, Q and z are cylindrical coordinates) 
and we have 



B^Bq Bg 



Tbrak 47rp/l Tbrak 



BzVe 



(1) 



where Tbrak is the braking time and h is the typical scale- 
height. Combining these two equations, we obtain 



^brak 



B, 



(2) 



which clearly shows that the braking time is shorter as, 
h, decreases. Assuming that the column density remains 
roughly constant, we find that Tbrak oc h^^^. 

As discussed in Galh & Shu (1993), Li & Shu (1996) and 
HF08, a magnetized cloud flattens along the field lines pro- 
ducing a thin magnetic pseudo-disk, which unlike centrifu- 
gally supported disks, is dynamically collapsing. We can 
relate the scale-height of the pseudo-disk, h, to the angle 
a. For an aligned system, assuming mechanical equilibrium 
along the z-axis, one can infer the relation (see HF08 for 
details) 



(3) 



which refiects that the gravity and magnetic pressure com- 
press the gas along the z-axis while the thermal pressure 
resists the compression. In this equation Br , and the gravi- 
tational force dz(j), are estimated just above the disk (z ~ h) 
while pc is the equatorial (z ~ 0) gas density. As shown in 
HF08, for values of /i lower than 5, the magnetic pressure 
produces a very significant density enhancement (typically 
a factor 10 or more) in the equatorial plane and very stiff 
gradients along the z direction. We now consider a cloud 
that rotates along an axis oz' , belonging to the x-z plane, 
such that the angle, a, between oz and oz' is small. We de- 
fine f2 the rotation speed. The rotation can be decomposed 
into two components; one is a rotation about the z-axis 
at a velocity cos(a) and the other about the x-axis at a 
(time dependent) velocity ~ f2sin(a). As a consequence, 
the magnetic field lines are stretched. As in the aligned 
case, a toroidal component is created by the rotation with 
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Fig. 2. \BzBe/ p\ (see Eq. [1]) as a function of z at ?/ = 300 
AU. Solid line is 6* = 0°, dotted line 6 = 10° and dashed 
line 6* = 20°. 

respect to oz but another component, By ~ B^^ sin(Q;)r^/, 
is also generated by the rotation with respect to ox, where 
Tff is the cloud freefall time. Therefore, unlike the aligned 
case, the components of the magnetic field parallel to the 
x-y plane do not vanish in the equatorial plane. Assuming 
that the angle a and the global cloud rotation are small, the 
magnetic field remains roughly unchanged with the excep- 
tion of its y-component. The mechanical equilibrium equa- 
tion becomes 

PcC^ + B^/STT-B^/STT-f (4) 

where Br is the radial field corresponding to the case with- 
out rotation. Equation 2] shows that the By component 
generated by the cloud rotation limits the magnetic com- 
pression along the z-axis. Therefore, the typical width of 
the pseudo disk, h, increases as a increases, and the effi- 
ciency of the magnetic braking (cx is consequently 
reduced. In addition, since the rotation axis is not parallel 
to the z-axis, the rotation itself will also tend to impede the 
contraction along the z-direction. 

3. Numerical simulations 

3.1. Numerical setup 

To carry out our numerical simulations, we used the 
AMR code Ramses (Teyssier 2002, Fromang et al. 2006). 
Throughout the simulations, the Jeans length is resolved 
with at least 10 cells. An HLLD solver is employed. The 
initial conditions consist of a 1 solar mass spherical cloud 
whose profile, given by p{r) = Pc/{^ + {fl^oY), resembles 
the observed cores. We impose a density contrast of 10 be- 
tween the central density = 8 x 10^ cm~'^ and the edge 
density, p^. Outside the cloud, the gas has a density of Pe/10 
and it is in pressure equilibrium with the cloud edge. As in 
HF08, we use a barotropic equation of state. The cloud is 
initially in solid body rotation and threaded by a magnetic 
field, along the z-axis, whose intensity is proportional to 
the total column density through the cloud. The rotation 
axis is in the x-z plane and at an angle a with respect to 
the magnetic field. The ratio of rotation to gravitational 
energy is equal to /? ~ 0.03, while the degree of magnetiza- 
tion is determined by the parameter p. The initial thermal 
to gravitational energy ratio is ~ 0.25. 

3.2. Results 

For various values of p and a, Fig. [1] shows the density and 
velocity distribution in the x-y plane (first row) and y-z 



plane (second row). These correspond to the initial direc- 
tion of the magnetic field being respectively perpendicular 
and parallel to the plane of the snapshot. Note that the 
rotation axis is also either perpendicular or parallel, except 
in the case of a = 20° for which it is slightly tilted. At the 
time of the snapshots, the mass of the gas denser than 10^° 
cm~'^ (the density at which the gas becomes adiabatic) is 
~ 0.15 M0, while the mass of the gas denser than 10^^ 
cm-3 is ^ 0.1 Mq. 

The first two columns show results for the cases p = 16 
and p — 5, respectively, with a = 0°. In the p = 16 case, 
a dense disk whose diameter is ~ 200 AU forms around 
the central object. The size of the disk is estimated by 
determining the position at which the velocity is roughly 
perpendicular to the disk radius, so that the rotation ve- 
locity dominates over the infall motions. The cut across the 
y-z plane also shows an outflow being launched perpendic- 
ularly to the disk. In the p — 5 case, no disk of signifi- 
cant size forms. The gas instead falls directly towards the 
central thermally supported core. The entire configuration 
is unstable leading to significant departures from the ax- 
isymmetric configuration. The cut across the y-z plane also 
shows the thin pseudo-disk (| y |> 200 AU, | z |~ 100 AU 
for p = 10^ cm^'^), which is compressed by the magnetic 
pressure exerted by the curved field lines. The top panel of 
the third column {p = 5, a = 20°) shows that a disk- like 
structure of ^ 150 AU forms around the central object (as 
shown by the velocity field, rotation dominates over infall 
in this structure). We note that the structure appears to 
be non-axisymmetric and that the entire velocity field is 
more complex than in the first case {p = 16, a = 0°). For 
p = 5 with a — 20° (bottom panel of third column) the 
pseudo disk, visible at | y |> 200 AU, is twisted and overall 
thicker by a factor of ~ 2 — 3 than the case p = 5, a = 0° 
(bottom panel of the second column) . Additionally, we have 
performed simulations of the case a = 10° for which we also 
find a disk. According to the analysis presented in sect. 2, 
we believe that the thickness of the pseudo-disk is at the 
origin of the formation of a rotationally dominated disk-like 
structure, since it reduces the magnetic braking efficiency as 
shown by Eq. ([1]). To confirm that the different behaviours 
are indeed related to differences in the magnetic braking 
efficiency, we plot in Fig. [2] the quantity \BzBg\/p which 
appears in Eq. The cases 6 = 0° (solid hue), 9 = 10° 
(dotted line) and = 20° (dashed line) are shown. As can 
be seen, the quantity \BzBe\/p is higher, and the two peaks 
are closer (equivalently the gradients between the peaks 
are stiffer), for the aligned case; implying that the quan- 
tity BzBq / [ph) is higher by at least ~ 50 — 100% than for 
the two misaligned configurations {Bg = at z = 0, for 
symmetry reasons) . This sudden change of behaviour as a 
increases, indicates, as suggested by Eqs. ([Sill]), that the 
aligned configuration is probably too restrictive for study- 
ing disk formation in collapsing magnetized cores. 

Figure [3] shows the density and velocity distributions in 
the y-z plane for larger values of a. The first panel illus- 
trates the p = 5 and a = 60° case, while the second panel 
shows results for a = 90°. In both cases a disk of ~ 150—200 
AU forms. We note that the densities in these disks are sig- 
nificantly lower than those observed in the ^ = 16 and 
a = 0° case (first row of Fig. [T]); their structure is clearly 
more complex, and we note the presence of extended spi- 
ral arms, which are almost one order of magnitude denser 
than in the rest of the disk. Finally, the third panel of Fig. [3] 
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Fig. 3. Density and velocity cut through the y-z plane. The centrifugally supported disk, when it exists, is seen nearly 
face on. First panel is for /i = 5 and a = 60°, second panel is for = 5 and a — 90° while third panel is for /i = 2 and 
a = 90°. 



shows results for a highly magnetized cloud rotating per- 
pendicularly to the magnetic field — 2, a = 90°). In 
this case, even though the magnetic braking is less efficient 
than in the aligned configuration, the magnetic field is so 
strong that no disk forms. Interestingly, we also see that a 
relatively strong outflow is launched. For completeness, we 
also considered the case {n = 3, a = 90°) for which we find 
that a small disk-like object of size 50-100 AU forms. Thus, 
we conclude that in the perpendicular case and within the 
framework of ideal MHD, disks can form for values of /i 
higher than ~3. It seems probable to us that, given the 
complex field geometry induced in these non aligned con- 
figurations, ambipolar diffusion as well as ohmic dissipation 
may play a role in the formation of small disks. 



4. Conclusion 

To study the formation of circumstellar disks, we have per- 
formed numerical simulations of the gravitational collapse 
of magnetized dense prestellar cores for a range of values 
of jj, and a (the initial angle between magnetic field and 
the rotation axis). Our main conclusions are the following. 
In the aligned case {a = 0°) disk formation is prevented 
for values of fi as high as ~ 5 — 10. In the non aligned 
cases, and even for relatively small values of a ~ 10 — 20°, 
disks can form for smaller values of /i ~ 4 — 5, correspond- 
ing to larger magnetic field intensities. The main reason for 
this change of behaviour is a decrease, for increasing a, of 
the magnetic braking efficiency, which is linked to an in- 
crease in the thickness of the magnetized pseudo-disk. This 
is caused by the magnetic field lines being twisted by the ro- 
tation about an axis parallel to the plane of the pseudo-disk. 
When a = 90° , we find that disks may form for smaller val- 
ues of /i, as long as /X > 2 — 3, and for even lower values of /x, 
disk formation does not seem to be possible. Nevertheless, 
for these highly magnetized configurations, the question of 
whether a disk may form at later times, or because of non- 
ideal MHD effects (Hosking & Whitworth 2004, Machida 
et al. 2008, Mellon & Li 2009), remains unanswered. We 
recall that although Belloche et al. (2002) observe a signif- 
icant amount of rotation in the envelope of IRAM04191, 
they exclude a disk of size larger than 20 AU. This result is 
very difficult to explain without a magnetic field, but agrees 
at least qualitatively with our conclusion if the source is 
sufficiently magnetized (Troland & Crutcher 2008). 



Another important point is that while in the purely hy- 
drodynamical case, and for the cloud parameters used, frag- 
mentation occurs; we do not identify fragmentation in any 
of the magnetized cases explored. It is however possible 
that with a better treatment of the thermal processes, the 
magnetized disks which form in the low magnetized config- 
urations, may also fragment. 
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Appendix A: Classical analysis of the magnetic 
braking: cloud embedded into an external 
medium 

In the context of interstellar clouds, the classical analysis of 
magnetic braking (e.g.; Mouschovias 1991, Shu et al. 1987), 
considers a rigid and dense axisymmetric cloud. We define 
Pc to be its density, R its radius and Z its height. The cloud 
is surrounded by a diffuse inter cloud medium of density pic. 
The typical timescale for magnetic braking Thr, corresponds 
to the time necessary for torsional Alfven waves induced by 
the twisting of the magnetic field line, to propagate over a 
distance I, such that the mass of gas swept by the waves is 
comparable to the mass of the cloud itself. At this point, 
a significant fraction of the cloud angular momentum has 
been transferred to the intercloud medium. 

Two cases can be considered. First, when the magnetic 
field and the rotation axis are aligned, the waves propagate 
along the magnetic field at the Alfven speed Va, and Ipic — 
Zpc leading to 



The second case corresponds to the magnetic field and the 
rotation axis being perpendicular to each other, and the 
waves propagating in the equatorial plane of the cloud. The 
intercloud medium, which at time t = Tbr is reached by the 
torsional Alfven waves, is located in a cylinder of radius I 
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and height Z. In this case, {(ji^ryaf' — ^i^^Pic — Pc, which 
gives 

^ ^J^ + l. (A.2) 

Va V Pic 

Since in typical astrophysical circumstances the inter- 
cloud medium has a density that is low with respect to the 
cloud density, Eqs. (lA.ip and (|A.2|) show that the braking is 
usually more efficient when the magnetic field is perpendic- 
ular to the rotation axis than when it is parallel. However, 
this conclusion is obviously correct only when i? ~ Z, i.e., if 
the cloud aspect ratio is not too different from 1 . In particu- 
lar, from Eqs. (|A.ip - (|A.2| ) we see that if Z/ R <C 1/ •\/ pd Pic, 
the braking time is shorter in the aligned case than for the 
perpendicular configuration. 

Previous studies have demonstrated that magnetized 
clouds are usually very flat because of the magnetic com- 
pression exerted by the radial component of the magnetic 
field (see e.g.; Li & Shu 1996, HF08). These magnetized 
sheets, which are called pseudo-disks, are perpendicular 
to the average magnetic field. In the same way, centrifu- 
gally supported disks are also very flat objects that are 
perpendicular to the rotation axis. For these extreme con- 
figurations, the magnetic braking time can obviously be 
longer when the magnetic field and the rotation axis are 
perpendicular than when they are parallel (depending on 
the respective values of i?, Z and pd Pic) - Note that strictly 
speaking, if the magnetic field and the rotation axis are 
not aligned with each other, the resulting structure is fully 
tridimensional rather than axisymmetric. 

It is interesting to compare Eqs. (jA.ip - (IA.2|) with 
Eq. (21). Even though the latter is identical to the first for 
Pc = Pic, there are two major differences. First, even when 
R c:i Z, the magnetic braking is not significantly more effi- 
cient in the perpendicular configuration than in the aligned 
one. The relative efficiency of the magnetic braking in the 
two configurations is, instead, directly proportional to the 
cloud aspect ratio. Second, as the cloud is compressed along 
the field lines, the quantity Z/^/p^ is simply proportional 
to Z^/^, implying that the magnetic braking time in the 
aligned configuration decreases. Again this is unlike the case 
of a rigid cloud embedded in a diffuse intercloud medium. 
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